Ratchetlike pulse controlling the Fermi deceleration and hyperacceleration 
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Using an ac driven asymmetric pulse we show how the Fermi acceleration (deceleration) can 
be controlled. A deformed sawtooth (Ratchetlike) pulse representing the moving wall in the static 
Fermi-Ulam model is considered. The time integral from the pulse over one period of oscillation must 
be negative to obtain deceleration and positive to obtain hyperacceleration. We show that while the 
decelerated case is chaotic, for the hyperaccelerated case the Lyapunov exponents converge to zero. 
Numerical simulations indicate that the hyperaccelerated case is ergodic in velocity space. Switching 
between different pulse deformations we are able to control the particle acceleration. Results should 
be valid for any pulse for which the time integral can be manipulated between positive and negative 
values. 
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Many years ago Enrico Fermi [l| suggested an acceleration mechanism of cosmic ray particles interacting with 
a time dependent magnetic field. Later on different versions of the original model for the Fermi acceleration were 
proposed, fn the first one, the Fermi-Ulam (FU) model, a bouncingparticle moves between a fixed surface and a 
parallel oscillating surface This model was shown to be chaotic In order to improve simulations a simplified 

version of the FU model was proposed called the static wall model. It ignores the displacement of the moving 
wall but keeps the essential information for the momentum transfer as the wall was oscillating. This static model 
was discussed in many aspects 0, 0, @, [3|i even for circular billiards Usually invariant curves in the phase space, 
found for higher velocities, prevent the particle to increase its kinetic energy without bounds. Recently the hopping 
wall approximation was proposed @, which takes into account the effect of the wall displacement, and allows the 
analytical estimation of the particle mean velocity. Compared to the simplified static model, the particle acceleration 



parti 

is enhanced. The second kind of Fermi accelerated model was proposed in 1977 by Pustyl'nikov [llj, who considered 
a particle on a periodically oscillating horizontal surface in the presence of a gravitational field. The above topics got 
attention in various areas of physics, ranging from nonlinear physics 0, i, MM 0, i, i [3l , atom optics [11, [11, [ll , 
plasma physics [l^, [3| to astrophysics [l7, 18, 19]. 

In this paper we use an ac driven asymmetric pulse to control the acceleration (deceleration) in the simplified FU 
model. The pulse is a deformed sawtooth driving law for the moving wall. This Ratchetlike pulse differs from the ac 
driven asymmetric pulses {symmetric sawtooth) used for the Fermi acceleration in the early work of Lichtenberg et 
al. [1] and proposed recently to control the motion of magnetic flux quanta [20| and to analyze the relative efficiency of 
mechanism leading to increased acceleration in the hopping wall approximation [loj . In the simplified Fermi model [1| 
the particle is free to move between the elastic impacts with the walls. Consider that the moving wall oscillates 
between two extrema with amplitude vq. The gravitational force is considered zero. The system is described by a 
map Mi(2)(V5i, 4>n) = {Vn+i, </>n+i) which gives, respectively, the velocity of the particle, and the phase of the moving 
wall, immediately after the particle suffers a collision with the wall. Considering dimensionless variables the FU map 
with the deformed sawtooth wall is written as 



Ml 



for (f>„ < Tji, and 



M2 



Vn+1 



Vn + —{(t>n - Vl) I 

m 

<Pn + M— FT mod [rji +772), 

Vn+l 



Vn + —{4>n - Vl) I 

m 

<Pn + fJ'—rp mod (771 +772), 

Vn+l 



(1) 



(2) 



for (jjn > 7^1, where n is the iteration number and /x is the maximum distance between the walls. Since in this simplified 
model the displacement of the moving wall is ignored, the modulus function is used to avoid errors due to successive 
collisions which may occur in the original model. In other words, if after a collision with the wall the particle continues 
to have a negative velocity (a successive collision will occur in the original model) , the particle moves beyond the wall. 
The modulus for the velocity injects the particle back and fixes the problem. 
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FIG. 1: (aj-fcj The shape of the pulses used in the simulations. The deformed sawtooth effect is obtained when Arj = 772—771 7^ 
0.0. For (a) A77 = we have the symmetric accelerated case A, (b) A77 < we have the deformed des accelerated case D and 
(c) Ari > the deformed hyperaccelerated case H. (d) Illustration of a mechanical device which mimics the effect of the soft 
version of the pulses from (b)-(c). The ratchetlike cylinder may rotate clockwise and counter-clockwise. 

wall (the pulse) for different values of the asymmetry A77 = 772 — The deformed sawtooth pulse is obtained 
when A77 7^ 0.0. Such pulses can be easily obtained from pulse generators. As we will see we obtain three kind of 
accelerations: For A77 = we have the symmetric accelerated case A [Fig. [H^a)]; For A77 < we have the deformed 
decelerated case D [Fig. [T]Jb)]; and for A77 > the deformed hyperaccelerated case H [Fig. [TJc)]. Figure [T]Jd) shows a 
very simple example of a mechanical device which generates a similar effect as the deformed sawtooth. Imagine many 
particles colliding elastically with the walls inside the square box, and a ratchet cylinder rotating at the center. The 
rotating ratchet cylinder mimics the wall oscillation of the deformed sawtooth from Fig. [T]^b)-(c). When an external 
motor starts to rotate the ratchet cylinder counter-clockwise, particles are hyperaccelerated. After that the motor 
starts to rotate clockwise and the h ype raccelerated particles are decelerated. A similar mechanical device was used 



for micromotors in a bacterial bath [21 1 



Numerically we reckon the average particles velocity at a given time n from 



1=0 ^ j=l 

where the index i refers to the ith iteration of the sample j, and ^ is the number of initial conditions. The average 
velocity (V) for the three cases A, H, D mentioned above are shown in Fig. [21 We iterate the map ([T|) or ^ for times 
n = 1 X 10^, and 3000 initial conditions in the interval < y < lO^'^ and < < (771 +772)- Figure [D shows the phase 
space (on the left) and mean velocity (on the right) for the cases A :A77 = 0.0 in Figs. [Ha)-(b); H -.Arj = 0.01 in 
Figs.[2lc)-(d); and D -.Arj = —0.01 in Figs. [2]Je)-(f). First observation is that the accelerated mode can be observed in 
the case A. As the particle velocity increases, regular island are observed and (V) increases slowly until ^ 10. These 
regular islands prevent the particle velocity to increase very fast. It is worth to mention that all initial conditions start 
inside the chaotic region at low velocities. The growth rate of (V) depends on the number of regular islands inside 
the phase space. This is an expected results since this case is similar to the simplified model studied by Lichtenberg 
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FIG. 2: Evolution of 50 chaotic orbits on the phase space V x (jj for the parameters /i — 10, vo = 0.2 and: (a) case A -.Ar] = 0.0, 
(c) H -.Ar] = 0.01 and (e) D :Ar] — —0.01. Mean values of the velocity (b), (d) and (f) calculated over 3000 trajectories related 
with the three case (a), (c) and (e), respectively. We call to attention that the order of magnitude from {V) changes drastically 
between cases A, H, D. 
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and Lieberman 0,01, where the harmonic force was considered. An interesting behavior appears when the deformed 
sawtooth is introduced. For a very small asymmetry A?/ = 0.01 (case H), the phase space is totally filled and no 
regular islands are observed [Fig.[2{c)]. The corresponding (V) increases very fast until (V) ~ 2 x 10^ [see Fig.md)], 
showing that the accelerated mode is enhanced (hyperaccelerated) when compared to the case A. Another kind of 
motion occurs when At^ = —0.01 (case D) where many regular islands appear in phase space [see Fig. HJe)] which 
again prevent the acceleration to increase without bounds, as in case A. The corresponding (V) remains almost 
constant for this example [Fig. [2l[f)]. A relevant property of case D is that: if the initial condition for the velocity is 
very high, the particle is decelerated. In other words, the hyperaccelerated particles from case D start to decelerate 
up to the time they reach the regular islands from Fig. [H^e), but now from above. Therefore, there is a regular island 
barrier which prevents particles to decelerated until zero for short times. For very long times, however, particles may 
pass between the regular islands and decelerate to zero velocities. Such motion close to regular islands may lead 
to "sticky" or tra ppe d [ 22l| trajectories which affect the velocity and the convergency of the finite time Lyapunov 
exponents [H, H S llj^ 

In order to explain the effect of the deformed sawtooth pulse, we rewrite maps ([1]) and ([2]) in the very high velocity 
regime (Vq 1). In this case the time evolution for the phase can be neglected and we have 



Vn+i « nFo 1 



Vim 



(4) 



From this expression for the velocity, we clearly see that when A?/ < the velocity is reduced by each iteration and 
the decelerated mode is obtained. When Ary = the velocity is constant and when Ar; > the velocity increases 
without bounds and we get the hyperaccelerated case. 

The above illustration shows the importance of a small deformation in the sawtooth wall motion on the dynamics in 
phase space and on the corresponding acceleration. Now we look at the stability of the trajectories for each situation 
A,H,D. Figure [3] shows the time evolution of the mean finite time Lyapunov exponent (LE), (A), calculated over 
1000 trajectories for the three distinct cases A, H, D. The LEs are determined using the Benettin's method 27[. The 
full line shows (A) for the case A, the long-dashed line for the case H and the dashed line for D. In the accelerate case 
A, (A) is constant since every initial condition (with low velocity) is trapped inside the chaotic region of phase space. 
For the case D, the (A) is almost constant, increasing very slowly. The time behavior of the mean LE in case D is 
similar to the mean LE from the A case in the following sense: the trajectory is trapped in the specific chaotic region 
in the phase space, below the islands of regularity [please compare Figs.^a) and (e)]. For the hyperaccelerated case 
H (long-dashed line) the mean LE decreases almost to zero, showing that the filled space phase observed in Fig. ^c) 
is not chaotic. Albeit being not chaotic, no regular islands are observed, as in the other cases. This is apparently 
a typical behavior of linear unstable but ergodic systems [i^l, where besides having zero LEs the velocity space is 
totally filled. 
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FIG. 3: Mean value of the Lyapunov exponent calculated over 1000 trajectories for the three situations: A : Arj = 0.0 (full 
line), D : A77 < 0.0 (dashed line) and H : Arj > 0.0 (long-dashed line). 

At next we show how the Ratchetlike time asymmetry can be used to control the Fermi acceleration in order to 
get a desirable velocity. Figure [H shows the time evolution of the particle velocity using a combination of cases H, D 
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FIG. 4: Velocity as a function of the iteration n. It shows how the particles velocity can be controlled by switching the sawtooth 
pulse appropriately between cases H (hyperaccelerated) and D (decelerated), and finally F (fixed wall). 



discussed above. While for the case H the hyperaccelerated mode is obtained, the case D is able to decelerated 
the particle velocity. When the desirable velocity is reached, the wall oscillation can be turned off (case F). By 
appropriately switching the Ratchetlike pulse between H, D modes, the velocity can be controlled. It is important 
to mention that the control is only possible due to the following particular property of the decelerated mode: the 
deceleration only occurs after the particle is hyperaccelerated, it does not occur for small velocities [see Fig. [^Teyr). 

To conclude, for a long time the Fermi acceleration has been studied in different models and applications [J, S, 
i, S, S, 0, S, S, M, El, Il3, US, 0, El, M, 113, 0, [3, li^. However, none of these works (but one 0) considered 



the possibility to control the particles in order to achieve a desirable velocity. Here we show the effect of a deformed 
sawtooth (Ratchetlike) pulse to control the particle velocity in the Fcrmi-Ulam model. Changing the asymmetry 
parameter from the Ratchetlike pulse we are able to get Fermi hyperacceleration and deceleration. In fact, the time 
integral of the pulse (the wall position) over one period of oscillation must be different from zero in order to get 
the control. When the mentioned integral is positive (negative), the Fermi hyperacceleration (deceleration) occurs. 
This allows us to, by switching the pulse between hyperacceleration and deceleration modes, to control the particles 
velocity. We also show that the decelerated and normal accelerated cases are chaotic, while for the hyperaccelerated 
case the LEs converge to zero. Additionally we give numerical evidences which indicate that the hyperaccelerated case 
is ergodic in velocity space. It would be interesting to observe the deformed sawtooth pulse to control the particle 
velocity in the presence of dissipation [2^ . 
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